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Abstract — Diversity analysis of a number of Multiple-Input 
Multiple-Output (MIMO) applications requires the calculation 
of the expectation of a function whose variables are the ordered 
multiple eigenvalues of a Wishart matrix. In order to carry out 
this calculation, we need the marginal pdf of an arbitrary subset 
of the ordered eigenvalues. In this letter, we derive an upper 
bound to the marginal pdf of the eigenvalues. The derivation 
is based on the multiple integration of the well-known joint 
pdf, which is very complicated due to the exponential factors 
of the joint pdf. We suggest an alternative function that provides 
simpler calculation of the multiple integration. As a result, 
the marginal pdf is shown to be bounded by a multivariate 
polynomial with a given degree. After a standard bounding 
procedure in a Pairwise Error Probability (PEP) analysis, by 
applying the marginal pdf to the calculation of the expectation, 
the diversity order for a number of MIMO systems can be 
obtained in a simple manner. Simulation results that support 
the analysis are presented. 

I. Introduction 

In wireless communications, the Wishart matrix arises from 
the MIMO transmission environment, where the channel ma- 
trix is modeled as complex Gaussian, as in the Rayleigh fading 
model [ 1 jj . In particular, if the channel matrix is available at the 
transmitter as well as at the receiver, the beamforming matrices 
can be obtained from the singular value decomposition (SVD) 
of the channel matrix to build the diagonalizing structure, 
known to be optimal to maximize the performance Q. In 
the uncoded version of this multiple beamforming scheme, 
the diversity order, an important performance measure of 
MIMO systems at the high signal-to-noise ratio regime, is 
determined by the subchannel with the smallest eigenvalue 
of the Wishart matrix Q, J4), 0. In general, the diversity 
order is calculated from the PEP, expressed as E where 
is a function of the I th ordered eigenvalue fii, and 
E[-] is the expectation operator. In order to carry out this 
calculation, we need to find the marginal pdf of the single 
eigenvalue /i; of the Wishart matrix. A first order polynomial 
expansion is used to derive the simple closed form expression 
of the marginal pdf in and (6), while the more accurate 
expression as the sums of terms of the form /ife 11 ^ 1 is provided 
in Q. The resulting diversity order of multiple beamforming 
with the I th eigenvalue involved is (TV— l+l)(M— l+l) where 
TV and M are the number of transmit and receive antennas, 
respectively [3|, [4|. 
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The average PEP between two codewords in the coded 
multiple beamforming scheme, on the other hand, requires the 
calculation of the expectation E [<fi(ni, ■ ■ ■ , fJry)], where <fi(-) 
is a function with multiple ordered eigenvalues involved Q. 
For the pairwise codewords whose corresponding function 
includes all of the singular values available from the SVD of 
the channel matrix, authors in [7 1 calculated the diversity order 
from the simple closed form expression of the average PEP, by 
making use of the fact that the sum of all ordered eigenvalues 
follows a chi-squared distribution. If </)(•) is composed of 
a subset of the ordered eigenvalues, the calculation of the 
expectation needs the marginal pdf of the eigenvalues. The 
closed form expressions of consecutive and an arbitrary subset 
of ordered eigenvalues are given in |8|, while the expressions 
for unordered eigenvalues are provided in (9] and fit)! . A 
difficulty exists in determining an analytical diversity figure 
with these prior approaches. They are typically in the form 
of a product of integrals to be calculated, and consist of the 
incomplete Gamma functions that enable numerical evaluation, 
but make the analysis difficult. 

In this letter, we propose a methodology to calculate an 
upper bound to the marginal pdf of the ordered eigenvalues. 
Then, we derive the diversity order by using the upper bound 
to the marginal pdf. Since the direct calculation of the marginal 
pdf is very complicated due to the multiple integration of the 
joint pdf which has the exponential function, we suggest an 
alternative function as a substitute for the joint pdf to simplify 
the multiple integration. The resulting diversity order is [N — 
Pi + 1)(M — pi + 1) where p\ is the index to indicate the best 
among the eigenvalues appearing in the </>(•) function. 
II. Problem Statement 

The elements of the MIMO channel H e £MxN are 
assumed to be Gaussian with zero mean and unit variance. 
In addition, the covariance matrix, which is defined as Vj = 
-E[hjhj] where hj is the j th column vector of H, and f 
stands for conjugate transpose, satisfies V 3 = I for all j. 
Based on the assumption above, the matrix HH^ is called 
uncorrelated central Wishart matrix |5|. The i th eigenvalue of 
HH', denoted by /i;, is sorted such that fii > \ij for i < j. 
Throughout this letter, we use X and Y as X = m&x(N, M), 
and Y = min(N,M). 

The average pairwise error probability that the receiver 
decides c instead of c as the transmitted signal is upper 
bounded by Q 



Pr (c ->• c) < E 



exp 



Y 



(1) 



2 



where 7 is the signal-to-noise ratio, and ctj is a given non- 
negative real value. We note that a bound of this form can 
be obtained for a number of MIMO SVD systems, e.g., ifTTl . 
H21, iTHl . |[T4l . Let's define a m i n as the minimum among 



the nonzero a values. Using the inequality J2j 
>J. ..„//.,. we rewrite ([TJ as 



> 



Pr(c ->■ c) < £ 



exp 




(2) 



where is the fc element of a vector p = [pi ■ ■ ■ pk] 
whose elements are the indices corresponding to non-zero a, 
i.e., a Pk ^ 0. Similarly, s = \s\ ■ ■ ■ s^ Y -k)] T is defined as a 
vector whose elements are the indices k such that a Sk = 0. The 
vectors p and s are sorted in increasing order. To calculate ||2), 
we need the marginal pdf of the K eigenvalues by calculating 
the multiple integration over the domain V s 



I Opi : 



P&i, ■■■ ,P-y) dp s 



dp Sl . (3) 



The joint pdf of the ordered strictly positive eigenvalues of the 
uncorrelated central Wishart matrices p ■ ■ ■ , py) in © is 
available in the literature (TJ, |[T5j as 

Y 

- Y. H 

p(pi,--- ,Hy) = ip(pi,--- ,Mi-)e 3=1 (4) 



where the polynomial ip (pi 



,p Y ) is 



^(Pi 



■^)=n^ Y n^-^) 2 



(5) 



Because we are interested in the exponent of 7, a constant 
multiplier, which appears in the literature, and is irrelevant to 
the exponent of 7, is ignored in ((5]) for brevity. 

The closed form expression of the marginal pdf can be 
calculated by evaluating (|3). However, this evaluation is 
complicated due to the multiple integration of the product 
of the polynomial and the exponential function in (O-©. 
Alternatively, we will now develop a method to get a simple 
expression for an upper bound to the marginal pdf. Then, we 
will use the upper bound to calculate (f2]). 

III. An Upper Bound to the Marginal PDF 

The complexity of the multiple integration to calculate 
the marginal pdf in (0) mainly comes from the fact that 
the elementary integration inside the multiple integration, 
fn y me ~ y dy, generates a large number of terms of the form 
x n e~ x for large m, i.e., 



y e 



<dy = [- 



y e »| 

-e~ x {x m + mx 
+ m(m — l)x' 



rn 
m— 1 
2 



(6) 



However, if we remove the exponential function from the 
elementary integration, the integration produces only one term, 
resulting in a much simpler multiple integration. In addition, 



since the eigenvalues of the Wishart matrix are positive and 
real, e~ Mi < 1 holds true for any i. This idea leads to a simple 
result of the elementary integration as 

1 



y m e~ y dy < 



~m+l 



(7) 

/o m + l 

To apply the idea above to the calculation of the marginal 
pdf, we introduce an alternative function 



,Py) 



if ai = (8) 



, ip(pi 



k = l 



, hy) e k=1 " if <y.\ > 

where the exponential factors irrelevant to the variables of 
integration are removed, except /ii which is kept in the case of 
a\ = 0. The reason for keeping /ii will be explained later. Cor- 
respondingly, let's define / (fi Pl , • • • , [i VK ) in a similar fashion 
to (O with p (/ii, • • • , py) replaced by p (pi, ■ ■ ■ , py), that 



is, 



/ {Ppi I ' ' ' I PpK ) 



pQii, ■ ■ ■ ,p Y )dp s 



dp Sl . (9) 



We see that p (pi, ■ ■ ■ , py) < p (pi, ■■■ , Py) for either the 
case of ai = or a\ > 0, and therefore, / (p Pl , • • • , p PK ) < 

/ (Mpi > • • ■ > Ppk ). where / (Ppi >■■■ 1 Ppk ) has a simpler ex- 
pression. We will employ © to calculate our bound in the 
next subsections. 

1) For a\ = 0: Since some factors of p (pi, ■ ■ ■ , py) are 
irrelevant to the variables of integration, they can be moved 
out. By defining a polynomial g (p Pl , p PK ) 

K K 

9 (M Pl , • • • , Vpk) = IJ P Pk T Y II (^p* -Vpj) , (10) 

fc=l jyk 

we rewrite (O as 



' Ppk ) — 9 (a*pi j 

^32-1 

x / e 



1 PPK ) e k 1 



-Ml 



I) 



(I 



,PY) 



9 (Mpi , ' ' ' , tipK ) 



dp s 



■dp S2 dpi (11) 



where f is a polynomial of the remained factors, 

9(Mpi,-" ,P-p k ) 

and dp Sl is replaced by dpi since s\ — 1 for a\ = 0. 



over the variables 



The multiple integration of —, v 

p S{Y K) , • ■ ■ , p S2 results in an intermediate polynomial whose 
terms are composed of the variables pi, p Pl , p PK . 
In other words, the intermediate polynomial is the sum of 
the terms of the form as p** x ••• x J Q p\e~^ 1 dp\. 
The final integration of the intermediate polynomial over pi 
leaves a polynomial with the variables p pi , ■ ■ • , p PK since 
L p\e~^ 1 dpx = y\. It should now be clear why we keep 
e~ Ml in p (fix, • • • , p Y ), The absence of this factor would have 
led the integration to diverge. 
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Defining hx(p pi , - • • , /Vk) as the result of the multiple 
integration, we rewrite ( fTTT i as 



1 I ' ' ' ) MpfC 



- E M Pfc 



(12) 



where the polynomial ri (/it Pl , • • • , ) is defined as 



T\ (Mpi i ' ' ' i t^pi< ) 

)• (13) 

2) For «i > 0: Using the polynomial g (/i Pl , • • • , [i PK ) in 
( [Tol l, we rewrite (O in this case as 

K 

f (A*Pi i ' ' ' i ^Pk ) = 9 (Mpi ! ' ' ■ j Mpk ) e 



/0 JO 

9 (a*pi i ' ' ' i Mpic) 



rf/i^r^, • • ■ dfj, S2 d(i sl (14) 



where f^ 1 '"' is the same kind of polynomial that is 

gut?,!,— ,t*p K ) 

described in the previous subsection. The multiple integration 
results in a polynomial of the variables fi Pl , ■ ■ ■ , fi PK without 
diverging since the last integration is not involved with infinity. 
By defining fi2 (fi pi , ■ ■ • , [i VK ) as the result of the multiple 
integration of ( TBI , we get / (// Pl , ■ ■ • , \i PK ) as 

K 

- E n Pk 

f(v P i,--- >Mpk) = r 2 (MpiT-' .MpjfJe fc=1 (15) 
where r 2 (/Lt Pl , ■ ■ ■ , fi PK ) is defined as 

T~2 (Mpi > " " " J MpK ) = 

1 Mpx 1 MPK )• (16) 

The polynomials n (/x Pl , • • • , /i PK ) and r 2 (/i Pl , • • • , \i PK ) 
are multivariate polynomials with many terms. It is worthwhile 
to focus on the smallest degree of the terms because it plays 
an important role in determining the behavior of (Q]) in the 
high signal-to-noise ratio regime. The mathematical analysis 
of the behavior will be described in Section [IV] with the help 
of the next Theorem. 

Theorem 1: The smallest degree of the multivariate poly- 
nomial n (/i Pl , ■ • • , [i PK ) or r 2 (/i Pl , • ■ ■ , n PK ) is (JV - pi + 
1)(M- Pl + 1)-K. 

Proof: See Appendix lAl □ 

In the case of the single eigenvalue where K = 1, and 
pi = I, Theorem Q] states that the smallest degree of fii is 
(N-1 + 1)(M - / + 1) - 1. This generalizes the result of the 
first order expansion in 0, (6) to calculate the marginal pdf 
of the I eigenvalue. 

IV. Calculation of the Expectation 

According to the analysis in Section [HI] the marginal pdf 
/ (fi pi , /J, PK ) is upper bounded by the general expression 



/ (A*Pi ) ' ' ' ! Mpir ) — r (Mpi : 



- E ^Pfc 

1 Mpjf ) e k 1 



where r (/x pi , • • • , fi PK ) is a polynomial with the smallest 
degree of (N ~ pi + 1)(M — p\ + 1) — if. We are now ready 
to obtain an upper bound to ® by calculating 




A*pif ) d[i PK ■ ■ ■ dfi Pl 

T (Mpi ! ' ' ' ) Mpjf ) 



(17) 



-(l+7Q mi „) J2 Pp k 

x e fc=1 d/i P/c • • • d[i Pl 

where T> p is the domain of integration. Note that 

e -(i+ 7 a ml „) £f =1 Mpfc < e - 7 a ro4 „ Ef =1 M Pfc _ In Theorem|2] we 

provide the result of the multiple integration of a term whose 
degree is Y.k Ph- 

Theorem 2: For a multivariate term with variables 0; for 
i = 1, • • • ,K whose exponent, denoted by /3j, is a non- 
negative integer, the multiple integration in the domain oo > 
0i > 2 > • • • > 6» A - > is 



^■■■9 K K e *=i dO K ---d6 l = 
i j i.i 

K 



V *=i ^ (18) 

where £ is a constant. 

Proof: See Appendix iBl □ 
Since the polynomial r (/i Pl , • • • , /i pjf ) is the sum of a 
number of terms with different degrees, the result of dTTb is 
also the sum of the terms of (ja m i n ) whose exponent obeys 
Theorem [2] For large 7, it is easy to see that the overall sum 
is dominated by the term with the smallest degree of r y~ 1 . 
Theorem [2] indicates that the smallest degree of j^ 1 results 
from the smallest degree of r (/i pi , ■ • • , [i PK ). Therefore, we 
conclude that (fl]i is upper bounded by 



< 



-(7V- P i + 1)(M- P i + 1) 



(19) 



where 77 is a constant, and irrelevant to 7. 

V. Simulation Results 

Fig. [T] shows the calculation of (Q]) where N = M = 3, with 
several specific a values through a Monte-Carlo simulation. 
The legend represents the a values as a vector notation 
[ai • • -as]- Three dotted straight lines are the asymptotes 
whose exponents correspond to 1, 4, and 9. The curves of 
the a values [1 0], [0.1 1], and [3 5], whose p\ is 1, are 
parallel to the asymptote of (3 - 1 + 1) (3 - 1 + 1) = 9. By 
comparing the slopes of the other curves with the asymptotes, 
we see that the analysis is supported by the simulation. 

Fig. [2] depicts the simulation result of N = M = 4. 
The dotted lines are the asymptotes of the exponents 4, 9 
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Fig. I . Simulation results together with asymptotic diversity orders that show 
the validity of the technique for N = M = 3. The vector [ai c*2 03] is 
provided in the legend. Note the importance of the first ct; equal to zero in 
determining the diversity order. 
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Fig. 2. Simulation results together with asymptotic diversity orders that 
show the validity of the technique for N = M = 4. Note the accuracy of 
the technique even if the first nonzero eigenvalue may be much smaller than 
the successive ones. 

and 16. A comparison of a slope with the asymptote reveals 
that simulation matches the analysis. Note that even though 
a corresponding to the smallest eigenvalue is 100 times that 
of the best eigenvalue, the slope is determined by the best 
eigenvalue. 

VI. Conclusion 

We derived an upper bound to the marginal pdf of the 
ordered eigenvalues of a complex central Wishart matrix. Such 
matrices arise in the analysis of MIMO SVD systems. Our 
bound employs an alternative function to simplify the multiple 
integration of the pdf. 



For MIMO systems employing SVD, a standard bounding 
technique provides a simple bound such as (QJ, see e.g., 
10, 13,(111, 021, O, Q31. By applying our result to the 
calculation of the expectation, one can calculate the diversity 
order (N — p% + 1)(M — p% + 1) in the high signal-to-noise 
ratio regime for a system where N and M are the number 
of transmit and receive antennas, respectively, and pi is the 
index to the smallest nonzero weight value in a PEP analysis. 
Simulation results provided here and elsewhere support the 
validity of the technique for a variety of MIMO SVD systems 
with a number of different parameters, see, e.g., 1111 . |12|, 

03, 01- 

Appendix A 
Proof of TheoremQ] 

1) For a\ — 0: Since v\ (n Pl , • • • jAW) is a product of 
two polynomials as shown in ( fT3b . the smallest degree of 
ri (n Pl , • • • , H PK ) is the sum of the smallest degrees of each 
polynomial g (fi Pl , • • • , \i PK ) and hi (fi Pl , n PK ). Let's 
define Dg tSma u est as the smallest degree of the polynomial 
g(fi Pl , ■ ■ ■ ,(ip K )- It is easily found that all of the terms in 
( TTOb have the same degree. Therefore, 

D a , sm aiiest = K(X -Y) + K(K - 1) (20) 

where the degree of K[X — Y) is contributed by the K factors 
of the form fj,p~ Y , and K{K— 1) comes from the ( 2 ) factors 

in the form of (fj, Pk — /i Pj ) 2 . 

To calculate the smallest degree of the polynomial 
hi (/i Pl , • • • , n PK ), we need to know the degree of the poly- 
nomial "pJiLLZ . The polynomial i/j (fj,i, • • ■ , fiy) in © 

has Y factors of the form [if~ Y and ( 2 ) factors of the 
form (/ij — fij) . The division by g (/j, pi , [i PK ) makes the 
common factors eliminated, leaving (Y — K) factors of the 
form nf ~ Y and ((^) — (^)) factors of the form (/i^ — /i^ ) 2 . 
Hence, the resulting polynomial f^ 1 '"' '^ y ) has degree 

D huorg = (Y — K)(X -Y) + Y(Y - 1) - K(K - 1) 

(21) 

for all of the terms of the polynomial. 

Obviously, there exists an integer e such that s e < pi < 
s € +i since p\ > 1. The integration over fij, for 1 < i < s e 
in ( flTT i makes these variables vanish because of the integra- 
tion to infinity due to fii, while that over the other fj,j for 
Se+i < j < S(y-k) converts those variables into the variables 
/j, Pl , • • • , fJ, PK ■ In the meanwhile, all the terms in f^ 1 '"' :Mr - > 

have different distributions on the degrees of the individual 
variables although they have the same degree as an entire 
term. Therefore, the smallest degree of hi (fj, Pl , • • • , [i PK ) is 
determined by the term which has the largest degree of those 
vanishing variables of pMliZ _ jt j s no i; necessary to find 

9{V P1 ,— -V-vk ) 

all the terms with the largest degree of the vanishing variables. 
Instead, we can see that one of those terms, whose degree is 
Dh 1 ,org, includes the factors 

i=l j>i 
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In this case, the degree corresponding to the vanishing vari- 
ables in ( 1221 is 

Dh 1: vanishing = (Pi ~ — Y) 

+ 2Y(p x -l)-p x (p x -l) (23) 

where (p x — 1)(X — Y) is contributed by the (p x — 1) factors 
of the form [if~ Y , and the rest of the degrees are calculated 
from the factors of the form /i?. 

Finally, the integration over fij for s e +i < j < sry-K) 
accumulates the degree of the current variables as well as 
the previous variables belonging to s. Hence, the degree of 
the variables /ifc for k > p\ is kept during the multiple 
integration. In addition, during each integration of \ij for 
•Se+i < j < s (y-k)7 the degree increases by 1 due to the 
fact that f^' [i? +1 d[ii+i = fi" + /(n + 1). Since s € variables 
from the original (Y — K) variables of integration vanished 
in h\ (/i Pl , • ■ ■ , li PK ), the degree to be added is 



D 



hi ,adde 



d = Y -K - Pl + l 



(24) 



where s € is replaced by (pi — 1). The smallest degree of 
r i (ftpi j * " * j ftpK ) can be calculated as 

Dri , smallest — ^g^smallest H - &h\, smallest 
— ^g^smallest 

H - \J-^hi,org & hi, vanishing ^hi, added) 



= {X-p x + l)(Y-p x + l)-K 



(25) 



where D 



hi,org 



-D 



hi ^vanishing 



stands for the degree which is 



kept during the integration over fij for s e + x < j < sty-K)- 
If X is equal to N, then Y is M, and vice versa, the smallest 
degree is (N - p x + 1)(M - p x + 1) - K. 

2) For a x > 0: Similarly to the case of a x — 0, the 
smallest degree of Ti (/j, Pi , [l Vk ) can be calculated in the 
same manner as in d25l l. The smallest degrees of (f20T > and 
(|2TT > apply to this case as well. However, since p x = 1 in 
this case, leading to p x < Si for i = 1, • • ■ , (V — K), no 
variable vanishes, resulting in Dh 2 , vanishing = 0. For the same 
reason, each of the (Y — K) variables of integration adds one 
degree. Therefore, Dh 2 , added = Y — K. By writing an equation 
similar to ( l25l l. we get the smallest degree of the polynomial 
r 2 j • ■ • i AW ) as 



D 



r-2i smallest 



= XY- K. 



(26) 



In general, the smallest degree of the polynomial 
r x (n Pl , [i PK ) or r 2 , • • • , p, PK ) can be expressed as 
(N —p x + 1)(M — p x + 1) — K since this holds true even for 
the case of a x > where p x = 1. 

Appendix B 
Proof of Theorem[2] 

The first integral for the variable 9k in ( fT~8l > can be calcu- 
lated as 



'K 



; d6 



K 



Pk+1 

E 



{p K -i + i)\ 



UJ 



'K-l 



(27) 



We will ignore all the constants for the simple expression since 
we are interested in the exponent of lu. The second integral 
can be calculated as 



K-2 rvK—1 



JO 



npK-1 n(>K -u(9 K -1+(>k) 



'K-l K 



d6 K d6 K „ 



EE 



a(PK+PK-i+2-i-j) -2u>8 K - 



'K-2 



'K-2 



(28) 



Even though the exact expression for the k th integral can 
be obtained by extending the procedure above, it is too 
complicated. A simpler method to calculate the exponent of 
lj^ 1 after the final integral is reached by observing the fact 
that the sum of the exponent of a; -1 and 0k-\ in ( |27| | is 
{Pk + 1) for an Y term. This fact can also be observed in d28b 
as (Pk + Pk-i + 2). By defining Tk for the k th integral in 
the same manner above, we can generalize Tk as 



Tfc = ffg-t+i + k. 



i=l 



The final integration is calculated by 

6^ f (w, 6 X ) e-" ei d6 x 



(29) 



(30) 



where 8 X ) is the (K — l) th integral of the variables uj and 
Q\ with r^-i- In other words, (f30b is the sum of the many 
terms which have the form 



UJ 



j(v+0O e -**fc d 8 x 



(31) 



where x + y = Yk-i, and z is a constant depending on the 
term. We can easily see that OTI ) results in 

5w -(*+*+ft+i) (32 ) 
where <5 is a constant, and the exponent of oj^ 1 is 



K 



x + y + p x + l = K + J2^- 



(33) 
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